Non-Gaussian Fluctuation and Non-Markovian Effect in the Nuclear Fusion Process: 
Langevin Dynamics Emerging from Quantum Molecular Dynamics Simulations 
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Macroscopic parameters as well as precise information on the random force characterizing the 
Langevin type description of the nuclear fusion process are extracted from the microscopic dynamics 
of individual nucleons by exploiting the numerical simulation of the improved quantum molecular 
dynamics. It turns out that the dissipation dynamics of the relative motion between two fusing nuclei 
is caused by a non-Gaussian distribution of the random force. We find that the friction coefficient 
as well as the time correlation function of the random force takes particularly large values in a 
region a little bit inside of the Coulomb barrier. A clear non-Markovian effect is observed in the 
time correlation function of the random force. It is further shown that an emergent dynamics of 
the fusion process can be described by the generalized Langevin equation with memory effects by 
appropriately incorporating the microscopic information of individual nucleons through the random 
force and its time correlation function. 
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The fusion of two nuclei is one of the major non- 
equilibrium processes in low energy nuclear reactions 
where the fluctuation and dissipation play important 
roles. It is rather difficult to describe the fusion pro- 
cess without significant simplifications. Under various 
assumptions, several macroscopic transport models have 
been introduced to evaluate the formation of a compound 
nucleus in heavy-ion fusion reactions [1, 2]. However, the 
microscopic mechanism on how two colliding nuclei fuse, 
especially how the relevant kinetic energy dissipates into 
the intrinsic degrees of freedom (DoF), remains a subject 
requiring further research. 

On the other hand, it is becoming feasible to get 
various microscopic information out of ab initio nu- 
merical simulations, like time-dependent Hartrce-Fock 
(TDHF) theories [3-7] , the many-body correlation trans- 
port (MBCT) theory [8], the quantum molecular dynam- 
ics (QMD) [9], the antisymmetrized molecular dynam- 
ics [10], and the fermion molecular dynamics [11]. The 
TDHF theory is mainly based on the mean-field concept; 
in TDHF, fluctuations of collective variables are consid- 
erably underestimated. Much effort has been made to 
give a beyond-mean-field description of fluctuations [12]. 
The n-body correlations are incorporated in the MBCT 
theory [8] which has only been used in very light sys- 
tems [13]. 

The QMD is a microscopic dynamical n-body the- 
ory which was successfully used in intermediate-energy 
heavy- ion collisions (HIC) [9]. An improved QMD 
(ImQMD) has been developed in order to extend the ap- 
plication of QMD to low-energy HICs near the Coulomb 
barrier [14]. A series of improvements were made in the 
ImQMD; in particular, by using the phase space occu- 



pation constraint method [15], the Pauli principle, which 
is very important for low-energy collisions, is effectively 
taken into account. Making full use of the microscopic in- 
formation provided by ImQMD simulations, in this Let- 
ter, we try to understand how the macroscopic fusion 
dynamics emerges out of the microscopic one. 

We focus on a simplest case of symmetric fusion pro- 
cess with the impact parameter equal to zero. In this 
case, the whole system can be divided into the left- and 
right-half parts instead of a projectile and a target be- 
cause many nucleons frequently transfer between these 
two nuclei after they touch each other [16]. The relative 
motion between two centers of mass (CoM) of the left 
and right parts is chosen to be the relevant DoF described 
by the Langevin equation. Our analysis is limited in a 
stage where the relative distance r is much larger than 
its width. 

The one-dimensional generalized Langevin equation 
with memory effects reads [17, 18] 



du(t) 
dt 



= - [ 7 (t - t')u(t')dt' + -R(t) - - 
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1 dV(r) 



dr 



(1) 



where u(t) is the relative velocity between the two parts, 
R(t) the random force felt by either part, and fi the re- 
duced mass of the system. 

In the ImQMD model [14], a trial wave function is re- 
stricted within a parameter space {r^, p^}, where Vj and 
Pj are mean values of position and momentum opera- 
tors of the jth nucleon which is expressed by a Gaussian 
wave packet. The time evolution of the trial wave func- 
tion under an effective potential is governed by the time- 
dependent variational principle [9-11]. An expectation 
value of the Hamiltonian is given by using an improved 
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Skyrmc potential energy density functional. In this Let- 
ter, we concentrate on head-on collisions of 90 Zr+ 90 Zr. 
Ten thousand collision events were simulated. Each sim- 
ulation is started at r = r = 30 fm and with an inci- 
dent energy E = 195 MeV which is a few MeV above 
the Coulomb barrier. Numerical details can be found in 
Refs. [19]. 

A potential for the relative motion is defined as, 

V(r) = E tot (r) - E left (r) - E Hght (r), (2) 

where -Etot, -EWt, and -E^ght represent the energy of the 
total system and those of the left and right parts, respec- 
tively, each of which consists of the kinetic energy, the 
nuclear and the Coulomb potential energies. This po- 
tential is shown in Figs. 1 and 3. The TDHF has also 
been used to extract microscopic interaction potentials 
between two nuclei [5, 20] which gave similar features as 
potentials from the QMD simulations presented here and 
in Refs. [21]. 

The random force in the ith event is defined as 

R(x)i=Fi(x)-{F(x)), x = torr, (3) 

where Fi(x) = Y^f=ifi{ x ) denotes the total force act- 
ing on the left (right) part of the whole system in the 
ith event, (F(x)) = ^ J27=i Fi( x ) the mean value of the 
force, and //(x) the force on the jth nucleon in the left 
(right) part. Here A means the number of nucleons con- 
tained in the left (right) part and n denotes the total 
number of events. In Eq. (3) and hereafter, (Q) denotes 
an average of Q over all events. 

Distributions of R(r) at several distances are shown in 
Fig. 1. The random force at r(t = 0) = ro shows a Gaus- 
sian distribution with the full width at half maximum 
(FWHM) r « 0.1 MeV/fm. In a region far away from the 
barrier, e.g., r ps 18 fm, the random force has a Gaussian 
distribution with T rj 0.5 MeV/fm. From a certain dis- 
tance, r ~ 13.5 fm, there appears a non-Gaussian shape, 
as is observed in Fig. 1. According to the shape of the 
distribution of R(r), one may divide the whole process 
into three regions. Region 1 represents an approaching 
phase up to the touching point: The distribution has a 
Gaussian form with rather stable and narrow width. Re- 
gion 2 is from the touching point to the barrier top: A 
non-Gaussian shape appears. Region 3 is from just inside 
the barrier top to the fusing phase: The distribution of 
R(r) has again a Gaussian shape with r w 15 MeV/fm 
which is almost two orders of magnitude larger than that 
in Region 1. 

To make clear what happens in Region 2, we divide 
the distribution of R(r) into a symmetric Gaussian and 
an asymmetric tail parts as is shown in Fig 2. We note 
that the Gaussian part has almost the same width as that 
in Region 1. The detailed structure of the random force 
can be studied by examining the strength and direction 
of the force felt by each nucleon. One typical event in 




Random force (10 3 GeV/fm) 

FIG. 1. (color online) Distributions of the random force R(r). 
Each inset shows the potential V(r) with the blue dot repre- 
senting the position where the system locates. The contour 
plots display the nucleon density distribution of the system. 



the symmetric part is shown in Fig. 2(a): All nucleons 
are well divided into two separated groups expressing the 
projectile and the target, respectively. Moreover, each 
nucleon locating in the left side of each nucleus feels a 
force toward the right (positive value), and that in the 
right side feels a force toward the left (negative value), so 
as to keep a stable mean- field. For all events in Region 1, 
the resultant force made by all nucleons in each nucleus 
is almost zero. Namely, the intrinsic structure of two 
fusing nuclei is kept almost unchanged, so is the width of 
the random force. This situation persists in events which 
belong to the symmetric Gaussian part in Region 2 as is 
shown in Fig. 2(a). 

A typical event in the asymmetric tail is shown in 
Fig. 2(b). Nucleons are roughly divided into two groups 
surrounded by solid lines. However, there appears a small 
third group within the dashed line. Since a few points in 
the negative (positive) force region express a set of nu- 
cleons which escape from the left (right) nucleus, and are 
being absorbed by the right (left) nucleus, a resultant 
force made by these nucleons gives a large right (lcft)- 
directed component to the random force. These trans- 
ferred nucleons move in an average potential formed by 
both the projectile and the target; they play a role to 
open a window. 

When the two nuclei come much closer, there occur 
more events which have more nucleons in the third group. 
Meanwhile, the other two groups, originating from the 
projectile and target, become closer to each other. Con- 
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FIG. 2. (color online) Distribution of the random force R(r) 
at r = 13.5 fm which is divided into the symmetric Gaus- 
sian (dark blue) and asymmetric tail (light blue) parts. Two 
typical events are shown in the inset: The abscissa and the 
ordinate express relative position Z of each nucleon and the 
force it feels in the Z direction. 
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FIG. 3. (color online) The correlation function (R(r)R(r)) 
[red dots, in (MeV/fm) 2 ] and the friction coefficient j(r) 
(blue squares, in 0.001 c/fm). The grey line shows the po- 
tential V(r). Pink diamonds represent (R(r)R(r)) calculated 
by eliminating events in the asymmetric tail. 



scqucntly, the asymmetric tail in the distribution of R(r) 
becomes larger. At the border between Regions 2 and 
3, it becomes very difficult to distinguish the event in 
the center part of the distribution from that in the tail 
part and all events are absorbed into a widely spreading 
Gaussian distribution. 

From the above discussion, it is concluded that the 
main microscopic origin of the random force, i.e., a two 
orders of magnitude enhancement of the random force 
is generated by individual nucleons in the third group. 
These nucleons also result in the abnormal behavior in 
the distribution of R(r), i.e., the long tail in Region 2 and 
a much larger width in Region 3 compared to Region 1. 

Next let us extract information for the macroscopic 
dynamics out of microscopic simulations. Assuming that 
the work done by the friction force is completely con- 
verted into the intrinsic energy Ei ntr (r), one gets the fric- 
tion parameter j(r) microscopically 

7W = « (4) 

(P)r 

withFf ric (r) = dE- mtr (r)/dr, E- mtT (r) = E tot (r)-E co n(r), 
and E co u(r) = p 2 /2fi + V{r). p denotes the relative mo- 
mentum between two CoMs and its mean value (p) r at a 
given r is defined as 

i 

where Pi(t) and ri(t) are the momentum and coordinate 
of the i-th event at time t and the following correspon- 
dence is used: For each event i, a time U is chosen in 
such a way that the relative distance takes a given value 
r, i.e., ri(ti) = r. By using thus defined tj, one may 



define r-dependent correlation function as 
1 " 

(R(r)R(r)) =-J2 Ri{U)Ri{U)\ {tAri(u)=r} ■ (6) 
»=i 

Although U needed to reach a given r is a little bit dif- 
ferent for different events, one may transform a function 
of the time t to that of the distance r for the ensemble 
of events by using this method. 

Figure 3 shows the correlation function (R(r)R(r)) and 
the friction coefficient j(r) which play decisive roles in 
the macroscopic description of dissipation phenomena. 
As is seen from Fig. 3, both of them take almost the same 
shape and their peaks locate at similar r. The friction 
coefficient of the fusion process induced by a head-on 
collision extracted from TDHF simulations shows similar 
strong peak structure just inside of the barrier [6]. 

To explore more deeply the dynamical relation between 
the microscopic motion of individual nucleons and the 
macroscopic dissipative motion, in Fig. 4 we show the 
time correlation function of the random force (R(r,r)) 
which is defined as, 

1 ™ 

(R(r,r)) = -^2MU)Ri(ti-T)\ {u[r(M=r} . (7) 
i=l 

Apparently, (R(v, t = 0)) = (R(r)R(r)). 

In Fig. 4 one clearly finds the non-Markovian effect. 
Especially when r = 12 ~ 10 fm, it is important to take 
account of memory effects generated by the microscopic 
motion of nucleons when one tries to properly evaluates 
macroscopic effects of the dissipation. Starting from the 
generalized Langevin equation (1) with memory effects, 
one gets the generalized fluctuation-dissipation (GFD) 
relation (R(0)R(t)) = fiksTj^t) which properly takes ac- 
count of the time correlation of the random force. There 
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FIG. 4. (color online) Time correlation function {R(r, r)} (7). FIG. 5. (color online) Effective temperature (8) and 



have been many ways to define the temperature for com- 
pound nuclei (see, e.g., Ref. [22]). Here we define an 
effective temperature for colliding systems by applying 
the GFD relation, 



non— Markov 



(r) 



<P>, 



1 



liks -Ffric(r) J 



dr(R(r,T)). (8) 



The effective temperature ^non-Markov obtained by 
taking an average over all events is shown in Fig. 5 by red 
dots. Here it should be noticed that T^on-Markov shows a 
peak around the range where the asymmetric tail appears 
in the distribution of R(r). Since the relative motion for 
events in the asymmetric tail is strongly affected by a 
few transferred nucleons between two fusing nuclei, i.e., 
by those in the third group of Fig. 2(b), one may deduce 
an important conclusion that the macroscopic dynam- 
ics of relative motion described by the one-dimensional 
Langevin equation (I) is not appropriate in Region 2. 
In other words, the appearance of the non-Gaussian dis- 
tributed random force indicates a necessity of introducing 
a new macroscopic DoF. Whether or not this new DoF 
may be related to the formation of a neck is an open 
question [1, 23]. 

After eliminating the events in the asymmetric tail in 
the distribution of R(r), one gets an effective tempera- 
ture I^orf-Markov wmcrl is depicted in Fig. 5 by pink dia- 
monds. The strength of R(r) after eliminating the asym- 
metric tail is also shown in Fig. 4. r^™_ Markov shows a 
consistent feature with \J Ei ntT which represents the tem- 
perature in the Fermi gas model. Note that according 
to our numerical simulations, the effective temperature 
obtained under the Markovian approximation is by an 
order of magnitude smaller than T non _Markov That is, 
the amount of energy lost from the relative motion and 
deposited into the intrinsic DoFs can be described by the 
generalized Langevin equation with memory effects, but 
not by a simple one under the Markovian approximation. 

In summary, we have discussed the generalized 
Langevin dynamics with memory effects by using both 
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the macroscopic and microscopic information extracted 
from ImQMD simulations for the fusion process. The 
central collision of the symmetric reaction 90 Zr+ 90 Zr 
with an incident energy just above the Coulomb barrier 
is studied in detail. It is found that the dissipation dy- 
namics of the relative motion between two fusing nuclei 
is associated with non-Gaussian distributions of the ran- 
dom force. In addition to the macroscopic information 
like the friction parameter 7(7") and the potential V(r) 
for the relative motion, the microscopic information of 
the random force as well as of its time correlation func- 
tion and a proper treatment of the non-Markovian (mem- 
ory) effect in the Langevin dynamics are decisive for the 
dynamics of emergence in the nuclear dissipative fusion 
motion. 
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